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ABSTRACT 

The pre-explosion observations of the type II-P supernovae 2006my, 2006ov 
and 2004et, are re-analysed. In the cases of supernovae 2006my and 2006ov we ar- 
gue that the published candidate progenitors are not coincident with their respective 
supernova sites in pre-explosion Hubble Space Telescope observations. We therefore 
derive upper luminosity and mass limits for the unseen progenitors of both these su- 
pernovae, assuming they are red supergiants: 2006my (logL/L Q = 4.51; m < 13M ) 
and 2006ov (log L/Lq = 4.29; m < 1OM ). In the case of supernova 2004et we 
show that the yellow super giant progenitor candidate, originally identified in Canada 
France Hawaii Telescope images, is still visible ~3 years post-explosion in obser- 
vations from the William Herschel Telescope. High-resolution Hubble Space Tele- 
scope and Gemini (North) adaptive optics late-time imagery reveal that this source 
is not a single yellow supergiant star, but rather is resolved into at least three dis- 
tinct sources. We report the discovery of the unresolved progenitor as an excess of 
flux in pre-explosion Isaac Newton Telescope i'-band imaging. Accounting for the 
late-time contribution of the supernova using published optical spectra, we calculate 
the progenitor photometry as the difference between the pre- and post-explosion, 
ground-based observations. We find the progenitor was most likely a late K to late 
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M-type supergiant of 8+j M . In all cases we conclude that future, high-resolution 
observations of the supernova sites will be required to confirm these results. 

Key words: stars: evolution - supernovae: general - supernovae: individual: SN 2004et, 
SN 2006my, SN 2006ov 



1 INTRODUCTION 

Type II-plateau (II-P) supernovae (SNe) have long been held to be the core-collapse induced ex- 
plosions of red supergiant stars - evolved massive stars with cool (~ 3000 K), expansive (~10 3 R©) 



Smartt 



20091) . That the progenitors of such 



hydrogen (H) rich envelopes (for a detailed review seeL 
SNe should be H-rich follows from the observatio n of strong H lines in the spectra of these events, 
the definition of a type II SN (see for example. lFilippenkolll997l) . For a star to undergo core- 
collapse it must be of sufficient initial mass in order to proceed to the ultimate stages of nuclear 
burning and thereby attain an iron core. Unable to counteract gravity through further nuclear fu- 
sion, the iron core contracts, but is supported against complete collapse by electron degeneracy 
pressure. When the core reaches the Chandrasekhar mass limit (~1.4M Q ), electron degeneracy 
pressure is overcome and the core rapidly collapses until neutron degeneracy pressure brings it to 
a halt, forming a neutron star. This collapse releases around 10 53 erg of gravitational potential en- 
ergy, which flows outwards from the collapsed core in the form of neutrinos. Roughly 1 percent of 
this neutrino flux (10 51 erg) couples with the outer layers of the s tar, causing t he supernova explo- 



sion. For a review of the core-collapse process see, for example, 



Bethe 



1990L Theory tells us that 



the initial mass threshold for a star to produce a collapsing core is between 7-12 M Q (|Heger et al. 



2003 



Eldridge & Tout 



2004; 



Siess 



20071 : 



Poelarends et al 



2008) 



Attempts to directly observe SN progenitors are hampered by a range of issues: the rarity of SN 
explosions (~1 SN per galaxy per 100 years); the difficulty in resolving single stars in galaxies 
outside of our own; and the lack of deep, high-quality pre-explosion observations of SN sites. Not 
surprisingly few SNe overcome these problems, but in recent years the archive of the Hubble Space 
Telescope (HST) has gone some way to increase the numbers that do. The diffraction limited view 
of HST has made it possible to identify single massive stars in galaxies at distances of up to ~ 20 
Mpc. 

For the last decade, several research groups around the world have been systematically search- 
ing for the progenitors of Core-Collapse SNe in archival imagery, for the most part in data from 
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HST (e.g.. 



Prieto et al. 



Van Dvk et al. 



2008 



2003a 



c, 



Botticella et al. 



Srnartt et al. 



2001 



Srnartt et al. 



2003 



Li et al. 



2005a. 



20071 : 



20091) . It was in archival HST and high quality Gemini tele- 



scope imagery that the first confirmed discovery of a red superg iant progenitor for a type II-P SN 



was made; that of SN 2003 gd, w ith an estimated mass of ~8 M Q (|Smartt et al. 



2004 



Van Dvk et al 



2003b 



Maund & Smart' 



itor of SN 2005cs (Maund et al. 



2009). This was followed by the detection of the red supergiant progen- 



2005 



Li et al. 



2006 



Eldridge et al.ll2007D . again in archival HST 



data. Several other SN II-P progenitor candid ates have bee n reported in the lit erature: SN 



(Mau nd & SmarttbOOa 



Van Dvk et a. 



2005al). SN 2006my and SN 2006ov (ILi et al 



( Elias-Rosa et al 



2003c ), SN2 004A fflendrv et al 



2005) 



20071). SN 2008 bk (IMattila et al 



2009J), and SN 2009kr (IFraser et al 



SN2004et(Lietal 



999ev 



2008), SN 2008cn 



20101) . All were identified in HST data, 



with the exception of the progenitor of SN 2008bk which was found in images from the VLT. 
More often than not, even where HST data were available, non-detections have been the norm. In 
these cases, qualified upper luminosity and mass limits have been set from the depth of the pre- 
expl osion images, the main caveat being the assumption that the progenitor was a red supergiant 



(e.g. 



Srnartt et al 



2003; 



Van Dvk et a 



In an extensive paper (|Smartt et al. 



2003c). 



OO9I) w 



itor population of type II-P SNe by examining the progenitor detections and detection limits of 20 
such events in a homogeneous fashion. From a maximum likelihood analysis of these mass limits 
assuming a Salpeter initial mass function (IMF), we found minimum and maximum initial masses 
for the progenitors of type II-P SNe of m min = 8.5+J 5 M Q and m max = 16.5 + 1.5M©. The minimum 
mass is in good agreement with that predicted by theory, but the maximum mass is lower than ex- 
pected. Stars with main-sequence masses of up to ~ 25 - 30M o are expect ed to retain most of t heir 



Srnartt et al. 



( 20091) we 



H-rich envelope, becoming red supergiants and exploding as SNe II-P. In 
estimate that, if our 20 progenitors were actually sampled from a population with masses between 
8.5 and 25M Q , we should have detected 4 progenitors of between 17-25M Q . That we detected none 
by chance has a probability of just 2 per cent. We have termed this 2.4<x result the red supergiant 



problem and refer the reader to 



Srnartt et al 



(|2009|) for detailed discussion. For several of our sam- 



ple of 20 objects in that paper we relied upon previously unpublished works by our group. In this 
paper we report our detailed analysis of the pre-explosion observations of three of these type II-P 
supernovae: SNe 2006my, 2006ov and 2004et. 



SN 2006my was discovered by K. Itagaki (jNakano & Itagak: 



20061) on 2006 Nov 8.82 UT (all 



times reported here are UT), located at ar J2 ooo = 12 h 43 m 40 s .74, 5 n m = + 16°23'14".l, some 27" 



west and 22.5" south of the centre of its host galaxy NGC 4651. 



Stanishev & Nielsen (2006) 



R. M. Crockett et al. 



spectroscopically class i fied SN 2006my as type II-P, simi 



(Elmha 



mdi et al. 



Li et al 



2003 



Hamuy et al. 



2001 



Leonard et al- 



ar to the class prototype SN 1999em 



20021 ) at 1-2 months past maximum. 



20071) have since shown that its discovery was probably ~3 months after explosion by 



comparing its light curve with that of SN 1999 em, and its spectra w ith with those of SN 2004dj. 



2006) on 2006 Nov 24.86 in M61 



SN 2006ov was discovered by K. Itagaki (|Nakano et al. 
(NGC 4303), located at a J2 ooo = 12 h 21 m 55 s .30, (5 J20 oo = +4°29'16".7, which is 5.5" east an d 51" 



Blondin et al. 



d2006l) clas 



north of its host galaxy centre. From a spectrum obtained on Nov 25.56, 
sified SN 2006ov as a type II, similar to the type II-P SN 2005cs around 1 month post-explosion 
although somewhat reddened in comparison. However, they do not specify the deg ree of red- 



Li et al. 



(2007) show that 



dening. Through detailed analysis of the SN light curve and spectra, 
SN 2006ov is not significantly reddened and was most probably discovered some 3 months post- 
explosion. Li et al. point out that by mistakenl y comparing its class ification spectrum with one of 



SN 2005cs at around 1 month post-explosion, 



Blondin et al. 



(|2006r) would have observed a much 



redder continuum for SN 2006ov and erroneously interpreted this as being due to extinction. Since 
type II-P SNe become progressively redder as they evolve, the older age at discovery could explain 

its colour. 

SN 2004et was discovered by S. Moretti dzwitter et al.l2004l) on 2004 Sep 27, located at a J20 oo = 20 h 35 m 25 s .4, 8 } 



in the nearby spiral galaxy NGC 6946. Spectroscopic and photometric monitoring revealed the 
presence of H Balmer lines in its optical spectra a nd a plateau lasting ~110d in its early-time 



optical lightcurves, confirm ing it as a type II-P SN (|Li et al 



2007 



Maguire et al 



Misra et al. 



2007: 



2005a; 



Sahu et al. 



2006 



Misra et al. 



201 Ol). Detection of t he SN at X-ray and radio wavelengths (|Argo et al 



2005 



Martf-Vidal et al. 



20071) suggested the presence of significant circumstellar ma- 



terial (CSM), most likely the slow, dense wind of the progenitor. Interaction of the fast-moving SN 
ejecta wit h this material would c reate a hot, shocked region producing radio synchrotron and X-ray 



emission (IChevalier et al 



2006). Furthermore, at about 3 years from the explosion the appearance 



of wide box shaped emission lines as well as rebrightening of the SN at mid-infrared wavelenths 
were interpreted as being due to the impact of the ejec ta on the progenitor CSM (Kotak et al. 



Zwitter et al 



(|2004r) . who measured the 



2009). Extinction towards SN 2004et was estimated by 
equivalent widths of the interstellar Na 1 D lines from a high resolution spectrum of SN 2004et 
taken shortly after disco very, and calculated a tota l reddening (host + Galactic) of B(B-V) = 0.41, 



using the calibrations of 



Munari & Zwitteil (11997b . It is to this value that the authors of the above 



refereed publications defer. 

Detections of candidate progenitors have been reported for all three of these type II-P events 
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SNe 2006my and 2006ov by 



Li et al 



Li et al. 



(2005a) 



(120071) (hereafter Li07), and SN 2004et by 
(hereafter Li05). In $2]and Owe argue that the candidate progenitors of SNe 2006my and 2006ov 
are not coincident with their respective SN sites. In §H1 our analysis shows that candidate, yellow- 
supergiant progenitor for SN 2004et is still visible some 3 years post-explosion. This yellow source 
is actually a blend of several point sources, which are resolved in new HST and Gemini adaptive 
optics images. We also present previously unpublished observations from the Isaac Newton Tele- 
scope (INT), in which we identify a red supergiant progenitor. 



2 THE PROGENITOR OF SN 2006MY 



A candidate red supergiant progenitor for SN 2006my was detected by Li07 in pre-explosion 
HST Wide Field and Planetary Camera 2 (WFPC2) observations (program GO-5375, PI: Rubin) 
of its host galaxy, NGC465 1 . They located the SN position on the HST data through alignment 
with ground based observations of the SN from the Canada-France-Hawaii Telescope (CFHT). 
Thirteen sources common to both a CFHT+MegaCam r'-band SN image and the WFPC2 F814W 
observation were identified and their positions used to calculate a transformation between the 
respective coordinate frames. The SN position was measured in the CFHT frame and transformed 
to the coordinate system of the archival HST data, yielding (x,y) coodinates of (410.61, 158.81) 
on WF2 of WFPC2 with a positional error of +0.45 pixels (+45 mas). This error was defined as 
the rms error of the transformation. Li07 identify a source within the error circle on the F814W 
image which they propose as the SN progenitor, characterising it as a single red supergiant star of 
main-sequence mass M = 10*fM©. 

We have since obtained HST+WFPC2 observations of SN 2006my as part of our HST program 
GO-10803 (PI: Smartt). The higher resolution of these data allow a more accurate alignment with 
the pre-explosion data to be determined. The details of our analysis are reported below. During 
our work, we became aware of a similar analysis usi ng our HST follow-up observations being per- 



formed contemporaneously by 



Leonard et al. 



(|2008r) . Both groups agreed that our analyses should 



proceed independently, and we have both arr i ved a t the same conclusions. Note that our basic 



result was originally published in 



Smartt et al. 



(120091) . which references this paper in preparation. 
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Figure 1. Aligned pre- and post-explosion HST+WFPC2 images of the field of SN 2006my. All the images are centred on the SN position (marked 
by the cross hairs) and oriented such that North is up and East is to the left, (a): pre-explosion F555W, (b): pre-explosion F814W and (c): post- 
explosion F814W with the SN clearly visible. A faint object is just visible close to the SN position in the pre-explosion F814W image. 



2.1 Astrometry and Photometry 

The pre- and post-explosion HST+WFPC2 observations of the site of SN 2006my were retrieved 
from the HST archivej at the Space Telescope Science Institute (STScI) via the on-the-fly re- 
calibration (OTFR) pipeline (see Tabled])- Each of the observations consisted of multiple expo- 
sures and these were co-added using the iraf stsdas task crrej in order to remove cosmic rays. 
The crrej task does not account for the significant geometric distortion suffered by WFPC2, 
which ranges from a few tenths of a pixel at the centre of each chip, to 2-3 pixels at the chip 



edges ( 



Gilmozzi et al 



1995 



Holtzman et al 



1995 



Casertano & Wiggs 



2001 



Anderson & King 



2003b . We have chosen to follow this reduction procedure (rather than employ the Multidrizzle 



image reconstructio n pipeline which corrects for geometric distortion (IFruchter & Hook 



Fruchter et al. 



20021 : 



2009)) in order to remain consistent with the reduction performed by Li07. In this 



way we are able to quote measured positions that can be directly compared with those of L i07. 
However, we go on to apply t he distortion corrections derived by Anderson & King ( 20031 ) and 



Kozhurina-Platais et al. 



(|2003f) before performing the alignment described in detail below. To cor- 
rect for the effects of geometric distortion on point source photometry, the crrej combined images 
were also multiplied by a geometric correction image. 

It was noted that the pre-explosion F555W and F814W images were well aligned, hence a trans- 
formation between the pre- and post-F814W images only was used to determine the position of 
the SN on the pre-explosion frames. The pre-explosion SN site was imaged on the WF2 chip of 
WFPC2 (pixel scale = 0.1"), while the SN was imaged on the PC chip (pixel scale = 0.0455") 
in our follow-up observations. Aided by the exquisite resolution of the HST data we identified 49 



http://archive.stsci.edu/hst/ 



2 flkl552bu - http://www.stsci.edu/hst/observatory/cdbs/SIfileInfo/WFPC2/WFPC2PixCorr 
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Table 1. Pre-and post-explosion observations of the site of SN 2()06my. 



Date 


Telescope+Instrument Filter 


Exp. Time 

(s) 




Pre-explosion images 




1994 May 20 


HST+WFPC2 F555W 


660 


1994 May 20 


HST+WFPC2 F814W 


660 




Post-explosion images 




2007 Apr 27 


HST+WFPC2 F450W 


1400 


2007 Apr 26 


HST+WFPC2 F555W 


1200 


2007 Apr 27 


HST+WFPC2 F814W 


1200 



common point sources between the pre- and post-explosion F814W images. The pixel positions 
were first correc ted for geometric d i stortio n usi ng the chip and wavelength dependent transforma 



tions detailed in 



Anderson & King! (120031) and 



Kozhurina-Platais et al. 



(2003). 



Anderson & King 



(|2003|) reports that these transformations are accurate to +0.01 pixels in the WF chips and +0.02 
pixels in the PC chip. The distortion corrected positions were subsequently used to derive a trans- 
formation between the "distortion-free" pre- and post-explosion coordinate frames. The transfor- 
mation (involving shifts, scales and rotations in x and y) was calculated using the iraf task geomap, 
and found to have a total (x and y rms combined) rms scatter of +0.19 WF pixels (19 masj^l 

To test the robustness of the fitted rms, we divided the 49 point sources into two lists, one of 
24 and one of 25 sources, and calculated two independent transformations between the pre- and 
post-explosion frames, finding rms errors of +0.19 WF pixels for each fit. Our aim was to test 
each transformation with point sources not used in its derivation. The transformed positions of 
such objects are not affected by over-fitting that may have occurred in the geomap calculation, and 
as such can be used to derive an independent and robust estimate of the transformation rms. 

To that end, we used the coordinate transformation derived from one list to map the post- 
explosion positions from the opposite list onto the pre-explosion coordinate frame. The fit rms 
was then calculated from the scatter of the mapped coordinates around the measured pre-explosion 
positions. In this way we measured independent, rms errors of +0.20 WF pixels for each of the 
transformations. These are consistent with the original rms values reported by geomap (+0.19 WF 
pixels), confirming there was little or no over-fitting (not surprising given that we employed only 
linear terms in our geomap fit) and that the geometric distortion correction was successful. We 
adopt the slightly larger +0.20 WF pixel (20 mas) total rms for all of the above transformations. 
(Details of all the transformation rms values are recorded in Table |2l) 



3 Unless otherwise stated, all coordinate uncertainties quoted it this paper are total values (i.e. x and y uncertainties combined in quadrature) 
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Table 2. Fit residuals for iraf geomap transformations between pre- and post-SN 2006my HST F814W observations. See text for details. 



49 source transform 
<r x = 0.124 
(7y = 0.148 
<T tou.1 = 0.193 

(residuals estimated from transformed positions 
of fit-independent source coordinates) 

<t x = 0.134 cr, = 0.140 

(Ty = 0.151 tr y = 0.145 

o-, ofo / = 0.202 <T, ota i = 0.202 



24 source tranform 
cr x = 0.125 
tr,. = 0.146 
ow = 0.192 



25 source transform 
cr v = 0.116 

£Ty =0.151 

ow = 0.190 



All values are in units of WF chip pixels (0. 1 " / pix) 



The position of SN 2006my in the post-explosion F814W frame was measured using the three 
centring algorithms of the daophot task phot - ce ntroid, Gaussian and ofilter - and the PSF-fitting 



pholQ (vi 



2000bl) . The mean of the four measure- 



ments was adopted as the SN position, SN(x,y) po .« = (416.08, 448.54) +0.05 PC pixels (2.5 mas), 



estimated from the range of the four measurements 



astrometric uncertainty of HSTphot ([Dolphin 



2000b). 



'his uncertainty is consistent with the limiting 



In order to map this SN position to the pre-explosion frame using our geomap transforma- 
tions, we first had to correct for geometric distortion. This wa s performed, as fo r the a lign- 



ment point sources, us i ng the distortion corrections detailed in lAnderson & King 



Kozhurina-Platais et al 



(2003) and 



(|2003r) appropriate for the PC chip and the F814W filter. The distortion 
corrected position was calculated as SN(x ,y ) post = (416.13, 448.68). 

This was subsequently transformed to the pre-explosion coordinate system using the 49 point 
source transformation, yielding a (distortion corrected) progenitor position of SN(x ,y ) / , re = (410.25, 
159.21) on WF2, with a total uncertainty of +0.20 WF pixels (20 mas). The transformation rms 
(20 mas) completely do minates over the contributions from the distortion correction (1 mas) 
(| Anderson & Kingl 120031) and the SN position measurement (2.5 mas). (We also applied the 24 
and 25 source transformations to the SN position, finding transformed coordinates within 0.01 
WF pix (1 mas) of our value of SN(x ,y ) pre above.) 

To compare our pre-explosion SN site position with that of Li07, we had to app 
ric distortion correction (WF2, F814W correction from lKozhurina-Platais et 
position, which was quoted for the distorted coordinate frame. The original Li07 SN site position 
(Li(x,y) pre = (410.61, 158.81) +0.45 WF pixels) becomes Li(x ,y ) pre = (410.62, 158.56). Com- 



y a geomet- 



20031) to the Li07 



4 Positions measured by HSTphot differ from those of daophot by -0.5 in both x and y coordinates. All positions quoted in this paper are corrected 
to daophot values 



paring with our value of SN(x ,y ) 
(75+49 mas). 



The progenitors ofSNe 2006my, 2006ov & 2004et 9 
= (410.25, 159.21) we find an offset of 0.75+0.49 WF pixels 



We remind the reader that the follow-up HST observations of SN 2006my were not available 
to Li07 at the time of their analysis, necessitating the use of ground-based CFHT data to perform 
the image alignment. The vast difference in resolution between HST and ground-based data makes 
it more difficult to identify common, isolated point sources that can be reliably used to align the 
images. This explains the larger Li07 transformation rms and the offset between our transformed 



positions. We therefore favour transformations performe d using the HS 



follo w-up data as being 



more accurate, both those presented in this paper and by [Leonard et al.l (|2008l) . This is in keeping 
with the work of all groups involved in SN progenitor searches, including Li07, who aim to obtain 
the highest resolution follow-up data possible to perform differential astrometry for precisely the 
reasons we state above. We stress that, had the SN 2006my HST data been available, we are 
confident Li07 would have obtained a similar result. 



Like Li07 we used the PSF-fitting photometry package HSTphot (version 1.1.7b) ((Dolphin 



2000bl) to perform photometry on the pre-explosion images. This software package includes sev- 
eral pre-processing tasks, which mask bad pixels, cosmic rays and hot pixels; determine the sky 
background and create co-added images from observations with the same pointing. Data input to 
these pre-processing tasks must be the original fits files as received from the HST archive. In order 
to compare our results directly with those of Li07 we were careful to choose exactly the same op- 
tions as detailed in their paper. During pre-processing we co-added the cosmic-ray split frames in 
each filter. When running HSTphot we chose Option 10, which turned on local sky determination 
and turned off aperture corrections since there were no good aperture stars. In such cases HSTphot 
applies default filter-dependent aperture corrections. We also matched our signal-to-noise (S/N) 
thresholds to those of Li07 - detection S/N threshold = 2.5cr; final photometry S/N threshold = 
3. Oct. 



Li07 detect a source of 5.6<x significance in F814W with a magnitude of m f814lv =24.47+0.2, and 
nothing in the F555W frame. We too find a source in F814W, with exactly the same magnitude and 
S/N ratio, which is indeed very close to our transformed SN position but is not coincident within 
the astrometric errors. (A source is also detected in F555W but is classified by HSTphot as having 
a flat, extended profile and is therefore rejected from the stellar photometry list). Having identi- 
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Figure 2. Magnified view of pre-explosion WFPC2 F814W image centred on the site of SN 2006my. The transformed SN position is marked by 
the circle as indicated, the radius of which corresponds to the \cr positional uncertainty. The position of the Li()7 candidate progenitor is marked 
by the smaller circle, the radius of which corresponds to the positional uncertainty we derive in this paper. This source is shown not be coincident 
with the SN. The image is oriented such that North is up and East is to the left. 



fled the source in F814W, we measured its position using the three centring algorithms^ within 
daophot (see above) and used these values along with the HSTphot measurement to calculate an 
average pixel position: Object(x,y) = (410.40, 158.70) +0.21 WF pixels ( 21 mas) estimated from 



the range of the four measurem ents. Correcting for geometric distortion (IKozhurina-Platais et al. 



2003 



Anderson & Kingil2003l) we found Object(x ,y ) = (410.23, 158.46). 



Comparing with our transformed SN position, SN(x ,y ) pre = (410.25, 159.21) + 0.20 WF pix- 
els, we find a separation of 0.75+0.29 WF pixels (75+29 mas)(see Figure [2]). This displacement 
suggests that the pre-explosion source and the SN are not coincident. 

We note that the original Li07 transformation, calculated using CFHT observations of SN 
2006my, did produce a SN position that was coincident with this source within the (larger) as- 
trometric uncertainties. Their identification of this source as a progenitor candidate was therefore 
perfectly valid at the time. For the reasons stated previously, we favour the differential astrome- 
try performed using the follow-up HST data as being more accurate. We therefore suggest that 
the Li07 candidate is most probably not the progenitor of SN 2006my given the di splacement we 



Leonard et al 



find be tween the SN site and the candidate source. This agrees with the findings of 
J2OO8). However, the ultimate test will be to re-image the explosion site after SN 2006my has 
faded to find out if the Li07 candidate, or indeed an y other object, has disap peared (similar to the 



analysis of the SN 2003gd progenitor carried out by 



Maund & Smartl 



2009). 



5 Note that the centring box used for the daophot centroiding task was 5 pixels wide, and hence the measured source position was not erroneously 
shifted towards the bright pixel to the south-east in Fig. [2] 
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2.2 Progenitor luminosity and mass limits 

Since we conclude that the progenitor of SN 2006my is most probably not detected in the pre- 
explosion observations, we must derive detection limits for these d ata in order to set lumin osity 



Maund & Smartt 



and mass limits on the unseen stellar pre-cursor (see, for example, 

signal-to-noise ratio, S/N, of a star imaged with a CCD is given by the equation 
p 

1 xtnr 



2005). The 



S/N = 



(i) 



+ <x 



bg 



where F star is the flux of the star in units of electron counts, e , and cr^ is the noise contribution 



from all other background sources other than the photon noise from the star itself ( ^F stur ). Here 
we define cr bg as 



«W = ^v,(l + £)(F«, + JP) 



(2) 



where n ap and n s u y are the number of pixels in the aperture and sky annulus respectively, R is 
the detector readout noise in e~, and F sky is the mean flux of the sky background per pixel, also in 
e~. The term (1 + n ap /n s ky) accounts for the noise incurred due to any error in the estimation of the 
sky background, F sky . This effect is mitigated by using large numbers of sky pixels to estimate the 
mean sky value. Rearranging eqn.Q]and solving the resultant quadratic gives 
(S/N) 



(S/N) + ^(S/N) 2 + 4<r 



(3) 



Taking values for all the input variables from the observations, we can use eqn. [3] to calculate 
the value of F star required to produce a detection of a desired S/N ratio. An apparent magnitude 
limit can then be calculated as 



m Um = -2.5 log 



+ ZP + apcor + CTE 



(4) 



K gain X exptime ) 

where gain, exptime, ZP, apcor and CTE are the detector gain, exposure time, zeropoint, aper- 
ture correction an d charge t ransfer efficiency correction respectively. The updated WFPC2 Vega- 
mag zeropoints of lpolphin ( 2000a ) are taken from Andrew Dolphin's website^. These are defined 
as the magnitude of a source producing a count rate of 1 Data Number per second (1 DNs 1 ), 



which is why we covert our flux to such units in 
mined from the tabulated encircled energy curves of 



quation [H Ap e rture corrections are deter- 



Holtzman et al. 



(| 19951) and CTE corrections 
(IDolphinll2000ah calculated using the latest equations from Andrew Dolphin's website 3 . 



6 http://purcell.as.arizona.edu/wfpc2_calib 
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Using an aperture of 2 pixels radius and a sky annulus of inner radius 10 pixels and width of 
5 pixels (all WF pixels), we calculated a 3cr detection limit for the progenitor of SN 2006my 
in the pre-explosion F&14W image: 3crFSl4W = 24.8. That we define a detection limit only in 
F&14W is due to the fact that we did not detect the SN progenitor. It is obviously impossible to 
obtain any colour constraints from a non-detection, and hence we must assume a spectral type 
in ord er to set luminosity and mass limits. In the case of a type II-P SN, such as SN 2006my 



(Li07; 



Maguire et al.ll2010h . we expect the progenitor to ha ve been a red supergiant star. S Ne II-P 



require progenitors of large radii of the order 400-1000Ro (|Chevaliei 



1976 



Arnett 



1980), which 



are typical for red supergiant stars. Of our pre-explosion observations the F814W data provides 
by far the most restrictive limits on such a red object. 



The apparent detection limit was converted to a bolometric magnitude using estimates of the 
distance and extinction tow ards the progenitor s tar, and bolometric and colour corrections appro- 



priate for a red supergiant. 



Solanes et al. 



(|2002|) have collected distance estimates for NGC4651 



from seven different sources, deriving a mean distance of d = 22.3 +2.6 Mpc (distance modulus 
jx = 31.74 + 0.25) and we adopt this value here. Li07 found no evidence for any host-galaxy ex- 
tinction in spectra of SN 2006my and hence, like them, we apply onl y a correction for the G alactic 



Cardelli et al. 



1989). Bolo- 



Drilling & Landoltl J2000) 



Maund & Smartd (|2005h 



extinction of E(B - V) = 0.027 assuming the reddening laws of 
metric and colour corrections for an M0 supergiant were taken from 
with a further correction between the F814W and Cousins / filters from 
An uncertainty of +0.3" 1 was assumed for the bolo metric correction, estim ated from the range of 
values for red supergiants of late K to late M-type (ILevesque et alj 120051) . The bolometric mag- 
nitude limit was found to be M bol = -6.13 + 0.39, which corresponds to a luminosity limit of 
log L/L Q = 4.35 + 0.16. Assuming the errors in the distance, extinction and bolometric correction 
are Gaussian in nature, one can set an 84 per cent confidence limit for the progenitor luminosity 
of log L/L @ =4.51. 



Comparing thi s limit with the final luminosities of Cambridge stars stellar evolutionary models 
Eldridge & Tout! (120041) shown in Figure [3] we can estimate an upper mass limit for the progenitor 
star. A detailed discussion of all aspects of our methodology, including its merits over compar- 
ing progenit o r phot ometry with the closest stellar model tracks on an HR diagram, is provided in 



Smartt et al 



(I2009() • Since metallicity affects the final luminosity of a star of given initial mass 



(see Figure [3b) we attempted to constrain the metallicity at the position of SN 2006my. The de- 
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Figure 3. Reproduction of Figure 1 from lSmartt et alj j2()09t) . (a): Initial mass versus final luminosity of STARS stellar models for solar metallicity 
(Z = 0.02). The luminosity range for each mass corresponding to the end of He-burning and the end of the model (just before core-collapse) are 
marked by the thin vertical lines. From a limit of luminosity an upper limit to the initial mass can be determined. The solid line is the luminosity 
of the model end-point, the dashed line the luminosity at the end of core helium burning and the dash-dotted line is the luminosity after second 
dredge-up when the lower mass stars become AGB stars, (b): The same as (a) but with three metallicities shown for comparison, and with the 
vertical joining bars omitted for clarity. 



UuJzl 



projected galactocentric radiusuof SN 2006my was found to be r G /r 2 5 = .37; almost exactly th e 



(2004). 



radius of the characteristic oxygen abundance measured in NGC465 1 by iPilyugin et al. 
They measure a value of log (O/H) + 12 = 8.7 dex at this position, which is consistent with so- 
lar metallicity (|Asplund et al.1 120041 8.66+0.05). We note that Li07 find an oxygen abundance 
of ~8.5 dex and co r rectly interpret this as subsolar compared to the solar value o f 8.8 dex from 



Grevesse & Sauval 



(|1998l) . With respect to the latest oxygen abundance estimates (|Asplund et al 



2004L 8.66+0.05) we interpret the Li07 measurement as being between our solar and LMC metal- 



licity values. From the LMC and solar metallicity models in Figure |3] we derive upper initial mass 
limits for the progenitor of SN 2006my of 1 1-13M Q . 



3 THE PROGENITOR OF SN 2006OV 

Li07 also reported the detection of a candidate red supergiant progenitor of main-sequence mass 
M = 15+3M for SN 2006ov in archival WFPC2 observations of M61 . In this case they aligned the 
pre-explosion fr ames with HST o bservations of the SN. Having performed PS F- fitting photometry 



using HSTphot (|Dolphinll2000b|) without detecting a progenitor star, Li07 noticed that a source 
was still visible in the residual images close to the SN site. They conclude that this residual object 
is coincident with the SN position, and report that by forcing HSTphot to fit a PSF at the SN 
position they detect an object of 6.1<x significance in the F814W and F450W observations. 

7 /'25 is the radius of a galaxy at which the surface brightness drops to 25 mag per square arcsec 
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Table 3. Pre-and post-explosion observations of the site of SN 2006ov 



Date 


Telescope+Instrument Filter 


Exp. Time 

(s) 




Pre-explosion images 




2001 Jul 26 


HST+WFPC2 F450W 


460 


1994 Jun 06 


HST+WFPC2 F606W 


160 


2001 Jul 26 


HST+WFPC2 F814W 


460 




Post-explosion images 




2006 Dec 12 


HST+ACS/HRC F435W 


420 


2006 Dec 12 


HST+ACS/HRC F625W 


180 
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Figure 4. Aligned pre- and post-explosion HST images of the site of SN 2006ov. All images are centred on the SN position (marked by the 
cross hairs) and oriented such that North is up and East is to the left, (a): pre-explosion WFPC2 F450W, (b): pre-explosion WFPC2 F606W (c): 
pre-explosion F814W and (d): post-explosion ACS/HRC F625W image with the SN clearly visible. The SN position lies just south-west of an 
extended source in the pre-explosion frames. 

3.1 Astrometry and Photometry 

We have analysed the same HST data-sets as detailed in Li07 and in Table [3] of this paper. The 
pre-explosion WFPC2 data consisted of F450W and F814W images from program GO-9042 (PI: 
Smartt) and F606W exposures from program GO-5446 (PI: Illingworth). The SN observations 
were taken with the Advanced Camera for Surveys (ACS) High Resolution Channel (HRC) (pixel 
scale = 0.025") as part of program GO-10877 (PI: Li). All data were downloaded from the HST 
archive at STScI via the OTFR pipeline, and the WFPC2 images further reduced using standard 
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tasks within iraf (see 32.11) . While the OTFR pipeline automatically corrected the ACS/HRC data 
for geometric distortion, the above reduction did not correct the WFPC2 data. Our data reduction 
was therefore entirely consistent with that of Li07, except in the case of the WFPC2 F606W 
pre-explosion observations. On inspection of this data-set we noticed that the pointings of the 
cosmic-ray (CR) split exposures were offset. With respect to the coordinate frame of the WF4 chip 
(on which the SN site was imaged), these offsets (exposure 1 minus exposure 2) were measured 
to be 6x = -0.12 pix and 6y = 1.16 pix. The CR-split frames were first aligned, using exposure 
1 as the reference frame, and then combined to produce the final reduced image (see Figure |4]). 
The "clipped" nature of objects along the y-axis of the co-added F606W image in Li07 (see Fig. 
8 in Li07), suggests that these offsets were missed. Li07 also used the HSTphot task coadd to 
combine the CR-split exposures before running the photometry package. This task cannot apply 
image offsets and will therefore have produced a similarly clipped image. Finally, the coordinate 
transformation we derive for this pre-explosion image results in a progenitor position that is offset 
in the y-direction from that of Li07 (see below). This is easily explained as a consequence of the 
differences in our respective reductions of the F606W data. 

We have repeated the alignment of the pre- and post-explosion HST observations originally per- 
fo rmed by Li07, except that here , as in 32. 1[ we have app l ied th e geometric distortion corrections 



of 



Anderson & Kingl (|2003|) and iKozhurina-Platais et al.l (|2003f) to the WFPC2 coordinates prior 



to calculating the transformation. 

In all of the pre-explosion observations the SN site was imaged on the WF4 chip of the WFPC2 
mosaic. The F450W and FS14W WFPC2 observations were found to be well registered and hence 
a single transformation was deemed appropriate for both. The positions of 3 1 sources common to 
both the pre-explosion F814W and the post-explosion F625W images were measured, corrected 
for distortion, and subsequently used to calculate a transformation (x and y shifts, scales and 
rotations) between the pre- and post-coordinate frames. This transformation had an rms error of 
+0.15 WF pixels (+15 mas). A separate transformation between the pre-explosion F606W image 
and the post-explosion frame was calculated from the positions of 31 common objects, with a total 
rms error of +0.15 WF pixels (+15 mas). 

We tested the robustness of the transformation rms values using the same method as employed 
in 32.11 The results of these tests are shown in Table HJ The residuals calculated from the fit- 
independent coordinates were only slightly larger than the geomap residuals, suggesting that there 



We assume a distortion correction appropriate for the F555W image 



Anderson & King 2003) 
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Table 4. Fit residuals for iraf geomap transformations between pre- and post-SN 2006ov HST observations. 



ACS/HRC F625W to WFPC2 F814W 

31 source transform 16 source tranform 15 source transform 

cr x = 0.081 cr x = 0.077 cr x = 0.085 

o\ = 0.122 (7y = 0.125 cr y = 0.106 

0"»m/ = 0.146 o-, ofo/ = 0.147 cr ma , = 0.136 

(residuals estimated from transformed positions 
of fit-independent source coordinates) 

cr x = 0.087 cr x = 0.079 

cr v = 0.137 <7y = 0.145 

O-wtal = 0. 162 CTtotal = 0. 165 



ACS/HRC F625W to WFPC2 F606W 

31 source transform 16 source tranform 15 source transform 

<t x = 0.095 cr x = 0.076 cr. r = 0.107 

cr y = 0.114 cr,, = 0.125 cr v = 0.100 

cr, OM ; = 0. 148 o-, om/ = 0. 146 a- tota i = 0. 146 

(residuals estimated from transformed positions 
of fit-independent source coordinates) 

cr x = 0.122 o- v = 0.091 

<r y = 0.107 cr v = 0.129 

0"»m/ = 0.162 tr Ioto/ = 0.158 



All values are in units of WF chip pixels (0.1"/ pix) 



was negligible over-fitting by the geomap task, and that distortion corrections had been successful. 
Nevertheless, we adopted the slightly larger +0.17 WF pixels (±17 mas) and +0.16 WF pixels (+16 
mas) as conservative, total rms errors for the WFPC2 F814W and WFPC2 F606W transformations 
respectively. 

The position of SN 2006ov in the post-explosion ACS/HRC F625W image was measured using 
the three centring algorithms within daophot and the PSF-fitting photometry package DOLPHOTOl 
The mean of the four measurements was adopted as the SN position, SN(x ,y ) poi; = (598.61, 
613.89) +0.12 ACS/HRC pixels (3 mas), estimated from the range of the four measurements. We 
again note that the ACS image was corrected for geometric distortion during pipeline reduction, 
hence no further correction was required. 

Applying the 31 source ACS-to-WFPC2 F814W transformation, we found a distortion corrected 
position for SN 2006ov on the pre-explosion FA50W and F814VK images of SN(x o ,y o )^ 50WV ™ 4H/ 
= (571.69, 236.19) +0.17 WF pixels (17 mas). This can be considered coincident with the trans- 
formed SN position from Li07, being offset by just 0. 12+0.24 WF pixels after correcting the Li07 
coordinates for geometric distortion (see Table[5]). 



9 http ://purcell. as . arizona.edu/dolphot/ 

10 Positions measured by DOLPHOT differ from those of daophot by -0.5 in both x and y coordinates. All positions quoted in this paper are 
corrected to daophot values. 
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Table 5. Comparison of transformed S N 2006ov positions i n pre- expl osion WFPC2 data from Li07 a nd this paper. Li07 coordinates corrected for 
distortion by applying the corrections of Anderson & King 12003) and Kozhurina-Platais et al. 12003). 



Pre-explosion 


Li07 position 


Li07 position 


Crockett et al. position 


image 




(distortion corrected) 


(distortion corrected) 


F450W/F814W 


(571.62, 236.22) ±0.17 


(571.76, 236.09) ±0.17 


(571.69, 236.19) ±0.17 


F606W 


(227.44, 267.34) ±0.16 


(227.64, 267.44) ±0.16 


(227.81, 268.01) ±0.16 



Coordinates are for WF4 chip (0.1" / pix), while uncertainties are total values (x and y errors combined) 



However, this is not the case for the F606W pre-explosion image. Applying the 3 1 source ACS- 
to-WFPC2 F606W transformation, we found the transformed distortion corrected SN position to 
be SN(x o ,y o )^ 06lv = (227.81, 268.01) +0.16 WF pixels (16 mas), which is offset by 0.59+0.23 WF 
pixels from the distortion corrected Li07 position (see Table [5]). This difference is due to the shift 
between the F606W CR-split images that was missed by Li07 during data reduction. We therefore 
adopt our F606W transformation as the most reliable. 

HSTphot was used to perform PS F- fitting photometry on the WFPC2 images. HSTphot and 
all pre-processing steps were run in precisely the same manner as in Li07 (see 92.11) . except in 
the case of the ^6061^ data where the offset between the CR-split exposures required special 
treatment. Here the exposures were not co-added; rather they were input separately to HSTphot 
along with appropriate offsets, using the first exposure as the reference frame. HSTphot returned 
no objects coincident with the SN position in any of the pre-explosion frames, but a point-like 
source was discovered close to the SN site in the FS14W residual image. Given that we have used 
identical HSTphot settings, we assume that this is the same F814W residual as observed by Li07, 
which they subsequently identified as a candidate progenitor for SN 2006ov. We, however, cannot 
conclude that this source is definitely coincident with the SN position. 

Figure [5] shows magnified sections of the WFPC2 pre-explosion frames along with the corre- 
sponding HSTphot residual images. All are centred on the transformed SN position, marked by a 
circle which represents the Icr astrometric uncertainty. 

The point-like F814W residual can be seen close to the centre of Figure [5]± We measured its 
position to be Res(x,y) f814W / = (571.16, 235.80) +0.12 WF pixels, which is marked in Figure [5]i 
by the smaller of the two circles just south of the SN position. Note that this is the raw position 
of the residual source, not corrected for geometric distortion. These coordinates can therefore 
be compared directly with those quoted in Li07. Comparing our residual position to the Li07 
transformed SN position, we found that the two are separated by 0.62+0.21 WF pixels. Correcting 



18 R. M. Crockett et al. 




Figure 5. (a-c): Zoom-in of SN 2006ov site on pre-explosion WFPC2 F8 14W, F450W and F606 W images. The transformed position of the SN is 
marked with a circle, the radius of which represents the ltr positional uncertainty, (d-f) Corresponding residual images, produced by HSTphot after 
subtraction of all detected sources. All images are oriented such that North is up and East is to the left. A point-like source is seen close to the SN 
site in the F814W residual frame (d). Its position, marked by the smaller circle, is not coincident with the SN within the combined uncertainties. An 
extended and less significant residual source is seen in the F450W frame (e). Again its centre is not coincident with the transformed SN position. 
No residual source is visible at, or near, the SN position in F606W (f). See text for details. 

the residual position for distortion, Res(x ,y )/r8i4w = (571.30, 235.67), and comparing with our 
own transformed SN position (Table© we found the same offset of 0.62+0.21 WF pixels. 

The F450W residual source (Figured) is neither point-like, nor as significant as that in FS 14 W. 
If it is a real source, its extended nature would suggest that it is not a single object. Nonetheless 
we have attempted to determine its centroid, Res(x,y)f 45(W = (571.17, 235.55) +0.13 WF pixels, 
which is marked on Figured and lies 0.81+0.21 pixels from the Li07 transformed SN posi- 
tion. Again, correcting the residual position for geometric distortion, Res(x o ,y o ) f 4 50 w = (570.86, 
235.91), and comparing with our transformed SN position we found an offset of 0.88 +0.21 WF 
pixels. 

The shallower F606W observation (Figures [Sfc&f) shows no apparent residual at, or near, the 
SN site. 

Li07 report that by forcing HSTphot to fit a PSF at their transformed SN position they de- 
tect a source of 6.1cr significance in F450W and F&14W, and 2.2cr in F606W. We performed 
the same forced-fitting, being careful to convert all daophot positions to HSTphot coordinates by 
subtracting 0.5 pixels from the x and y values. We were unable to reproduce their photometry at 
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the transformed SN site, finding detections of 3.6<x and 2.9a significance in F%\AW and FA50W 
respectively, and nothing in F606W. 

Detections of the highest significance (F450W = 24.16+0.24 (4.5cr); F814W = 23.36+0.18 
(6. Oct)) were found when forcing PSF fitting at our Residual source positions, which are offset 
from the transformed SN position. Admittedly, it is difficult to judge by eye which of the two 
circles, the SN position or the Residual centroid, is closer to the true centre of the residual source 
in the panels of Figure|5] However, since our best-fit photometry arises when fitting at the measured 
Residual coordinates, we favour this as its true position. 

Since we cannot reproduce the photometry of Li07 at the transformed SN site, and have found 
that the residual source is offset from the SN position, we suggest that this object is probably not 
the progenitor star. Given its proximity to the SN site, it is likely that some of the pre-explosion 
source flux in contributed by the progenitor. However, the evidence shown here suggests that it 
would be wrong to attribute all of this flux to the progenitor. As the residual itself is close to 
the noise limit in both FA50W and FS14W, we suggest that the progenitor of SN 2006ov is not 
observed. The ultimate test will be to re-image the explosion site with HST after SN 2006ov has 
faded to find out if the Li07 candidate, or indeed any other object, has disappeared. 



3.2 Progenitor luminosity and mass limits 



Luminosity and mass limits were set for the assumed red supergiant progenitor following the 
method described in §12.21 The sky background at the SN site was estimated and equations [2l|3]& 
|4]used to calculate the 3cr limiting magnitude of 3crF814V7 = 24.2. 

Estimates of the distance, extinction and bolometric and colour corrections were applied to 
convert this detection limit to a bolometric magnitude. Li07 calculated a mean distance m odulus 
for M61 of ix = 30.5 + OA (d = 1 2.6 + 2.4 Mpc) from two Tully-Fisher distance estimates (|Tully 



1988 



Schoeniger & Sofue 



1997|) . They also found no evidence for any host-galaxy extinction 



in spectra of SN 2006ov, applying only a small correction of E(B - V) = 0.022 for Galactic 
extinction. We adopt both these values here, and make the same assumptions for the bolometric 
and colour corrections as we have done in §12.21 In this way we derived a 3<x bolometric magnitude 
limit of M bo i = -5.48 + 0.50, corresponding to a luminosity limit of log L/L Q = 4.09 + 0.20. 
Assuming the errors are Gaussian, an upper luminosity limit of 84 per cent confidence can be set 
of log L/L Q = 4.29. 



The oxygen abundance gradient of M61 was re-defined by lPilyugin et al. 



(|2004l) and we use this, 
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along with the galactocentric radius r G /r 2 5 = 0.26, to estimate an oxygen abundance of 8.9 dex at 
the position of SN 2006ov. Hence, we compare the above luminosity limit with stars models of 
solar metallicity in Figure [31 finding an upper initial mass limit for the SN progenitor of 10M Q . 



4 THE PROGENITOR OF SN 2004ET 

A candidate progenitor star for the type II-P SN 2004et was identified by Li05 on archival images 
from the Canada France Hawaii Telescope (CFHT). These included BVR CFH12K images from 
2002 and u'g'r' MegaCam observations from 2003. The pre-explosion images were aligned with 
SN observations from the Katzman Automatic Imaging Telescope (KAIT), yielding a transfor- 
mation between the pre- and post-explosion coordinate systems with a rms error of ~ 0.1". The 
candidate progenitor was found to be coincident with the SN position to within these astrometric 
errors. Through comparision with stellar evolutionary models Li05 characterised this source as a 
single, yellow supergiant star with an initial main sequence mass of 15+5M©. Both the mass and 
the colour of this candidate progenitor were of note; the mass since it was the largest detected 
for a normal II-P SN, and the colour since such progenitors were expected to be red not yellow 
supergiants. 



To m ore accurately determine the position of SN 2004et on the pre-explosion frames 



Li et al. 



(I2005bh aligned the CFHT data with HST ACS/HRC observations of the SN. They reported that 
the transformed position of the SN was "consistent with the progenitor in the CFHT image to 
within 0.056" + 0.043"" and that this confirmed their identification of the progenitor star. 

We have since obtained deep follow-up observations of the site of SN 2004et from the William 
Herschel Telescope (WHT), the HST and the Gemini (North) observatory, which show that the 
Li05 candidate progenitor is still visible approximately 3 years post-explosion, and that it is in 
fact resolved into at least 3 separate sources. We also present hitherto unpublished pre-explosion 
/'-band observations of the SN site, taken at the Isaac Newton Telescope (INT), in which we detect 
an excess flux that we attribute to the unresolved progenitor of SN 2004et. 



4.1 Data acquisition and reduction 

The pre- and post-explosion observations for SN 2004et are detailed in Tabled The Canada France 
Hawaii Telescope (CFHT) datasets were downloaded from the CFHT archive at the Canadian 
Astronomy Data Centre (CADCHJ along with appropriate flat field and bias frames. The B&V data 

1 ' http://www4.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc 
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Table 6. Pre-and post-explosion observations of the site of SN 2004et. 



Date 


Telescope+Instrument 


Filter 


Exp. Time 

(s) 


Pre-explosion images 






2002 Aug 06 


CFHT+CFH12K 


g 


450 


2002 Aug 06 


CFHT+CFH12K 


v 


300 


2000 Sep 27 


CFHT+CFH12K 




360 


2002 Aug 12 


INT+WFC 




3600 


Post-explosion images showing bright 






2005 May 02 


HST+ACS/HRC 


F625W 


360 


Late-time post-explosion images 






2007 Aug 12 


WHT+AUX 


g 


2400 


2007 Aug 12 


WHT+AUX 


v 


1800 


2007 Aug 12 


WHT+AUX 


R 


1800 


2007 Aug 12 


WHT+AUX 


I 


2400 


2007 Jul 08 


HST+WFPC2 


F606W 


1600 


2007 Jul 08 


HST+WFPC2 


F814W 


1600 


2007 Jul 08 


HST+NICMOS 


F110W 


640 


2007 Jul 08 


HST+NICMOS 


F160W 


512 


2007 Jul 08 


HST+NICMOS 


F205W 


576 


2008 Jan 19 


HST+WFPC2 


F606W 


1600 


2008 Jan 19 


HST+WFPC2 


F814W 


1600 


2008 Jan 19 


HST+NCMOS 


F110W 


640 


2008 Jan 19 


HST+NICMOS 


F160W 


512 


2008 Jan 19 


HST+NICMOS 


F205W 


576 


2008 Jul 10 


Gemini-N+Altair/NIRI 


K 


7200 



were the same as those analysed by Li05. However, the R -band image was different. Taken in 2000, 
the total exposure time of 360 sec is only slightly longer than that of the 2002 7?-band image (300 
sec) analysed by Li05, but the image quality is 0.6" (compared to 0.8") and the SN site happened to 
fall on one of the high resistivity (HiRho) bulk silicon chips of the CFH12K mosaic. The sensitivity 
of these chips in the 7?-band is on average 15 percent higher than the standard epitaxial silicon (EPI) 
chips that make up the rest of the mosaic, and on which the observations from Li05 were taken. 
The net result is a deeper pre-explosion image than has yet been published. All CFHT data were 
reduced using standard tasks within iraf. In addition we discovered a deep INT+WFC i'-band 
observation from 2002, not analysed by Li05, which extended the pre-explosion SED to longer 
wavelengths. The INT data were downloaded from the CASU INT Wide Field Survey (WFS)l3 
having been processed and calibrated via the INT WFC pipeline. The HST Advanced Camera for 
Surveys/High Resolution Channel (ACS/HRC) observation of SN 2004et from 2005 was retrieved 
from the HST archive at STScI via the OTFR pipeline. The pipeline Multidrizzle task combined 
the sub-exposures and corrected the resultant image for geometric distortion. 
We re-observed the field of SN 2004et using the Auxilary Port Imager (AUX) on the WHT in 



http://www.ast.cam.ac.uk/ wfcsur 
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Aug 2007, ~3 yrs post-explosion. By this late-epoch we judged the SN would have faded suffi- 
ciently to allow us to check whether or not the Li05 candidate progenitor had disappeared; a key 
test of its validity. The observations were taken in BVRI filters to match as closely as possible the 
pre-explosion CFHT and INT images. All WHT data were reduced using standard tasks within 

IRAF. 

In an attempt to resolve the explosion site at late-time, after the SN had faded, ground-based 
adaptive optics (AO) observations of the field of SN 2004et were taken on 2008 July 10 (~ 4 yrs 
post-explosion) using Altair/NIRI on the 8.1-m Gemini (North) Telescope (programs GN-2008A- 
Q-28, PI: Crockett and GN-2008A-Q-26, PI:Gal-Yam). It was only afterwards that the authors 
discovered proprietary HST WFPC2 and NICMOS data, taken 2007 July 8, had become public 2 
days earlier. Further HST data, taken 2008 Jan 19, became public in Jan 2009. 

The Gemini observations were taken using a nearby 2MASS source (2MASS 20352132+6007102) 
as the natural guide star (NGS) for the Altair AO wave front sensor. This guide star was separated 
by around 30" from the SN position. Due to the geometry of the detector and wavefront sensors, 
the maximum displacement of a guide star from the centre of the FOV is 25". The site of SN 2004et 
was therefore imaged off-centre using the f/14 camera, which has a FOV of 51" x 51" (pixel scale 
of 0.05"). The effects of anisoplanatism, which would seriously reduce the image quality at such 
large separations from the guide star, were mitigated by inclusion of the field lens in the optical 
path. A series of short exposures totalling 7200 sec on-source were taken in the Z-band and sub- 
sequently reduced, sky subtracted and combined using the NIRI reduction tools within the iraf 
gemini package. The quality of the reduced image at our target position was 0.15"; excellent given 
the large displacement from the AO guide star and that the guide star itself was relatively faint for 
its purpose. 

The HST WFPC2 and NICMOS observations were carried out as part of program GO- 11 229 
(PI: Meixner). All data were downloaded from the STScI archive via the OTFR pipeline. WFPC2 
observations were taken in F606W and F814W filters, in each case as a series of 4x400 sec expo- 
sures. NICMOS Fl 10W, F160W and F205W data were taken as 5x128 sec, 4x128 sec and 4x144 
sec exposures respectively. Both the WFPC2 and N ICMOS sub-exposures were dithered, and we 
used the drizzle technique (|Fruchter & Hookll2002|) to combine the sub-exposures in each filter. 
This reduction process also corrected the WFPC2 and NICMOS images for geometric distortion. 
The SN site was imaged on the PC chip of the WFPC2 instrument (drizzled pixel scale = 0.0455") 
and on the NIC2 camera of NICMOS (drizzled pixel scale = 0.075"). 

In the following subsections we discuss first our natural- seeing, ground-based observations, 
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following by the HST and Gemini images, before drawing final conclusions from the combined 
data set. 



4.2 Ground-based, seeing-limited data: image alignment and photometry 

In Figure [6]we show the CFHT+CFH12K B&V images, as presented by Li05; the previously un- 
published CFHT+CFH12Ki? observation from 2000; and the deep INT+WFC i'-band observation 
from 2002, also previously unpublished. In the same figure we present our late time WHT+AUX 
BVRI observations, with the BVI images combined in Figure [7J to create a pseudo-colour image. 
It is immediately apparent that an object is still visible at the position of the Li05 candidate pro- 
genitor in all bands. Alignment of the CFHT and WHT R -band images demonstrated that the pre- 
and post-explosion sources are coincident, having a measured separation of 5 + 50 mas. 

PSF-fitting photometry of the Li05 candidate progenitor source was carried out on both the 
archival and late-time observations using tasks within the iraf daophot package. Although pho- 
tometry for the CFH12K B&V images was presented by Li05, it was decided to repeat this analysis 
so as to obtain measurements in a consistent fashion across all our pre- and post-SN observations. 
Zeropoints and colo ur corrections were derived using stars in the field with calibrated magnitudes 



taken from Li05 andlMisra et al. 



(120071) (see Figure[8]and Table©. The INT+WFC i'-band obser- 
vation was an exception, the zeropoint being taken from the image header (calculated during the 
WFS calibration process) and colour transformations to the Landolt system from the INT Wide 
Field Survey (WFS) websitef]. A series of Landolt standard stars were observed every night during 
the INT WFS, and secondary standards were also established in each field. The photometric cali- 
bration of the WFS data is therefore very reli able. Nevertheless we have checked this calibration 



using the secondary standards from Li05 and iMisra et al.1 (|2007l) presented in Table [/J and found 



that the results are consistent. Several of the isolated standard stars (Figure [8]) were used to build 
emperical PSFs for each image. Sources close to Li05 candidate progenitor were fitted simultane- 
ously, but their photometry is not detailed in this paper. Our pre- and post-explosion photometry 
is recorded in Table [81 along with the pre-explosion photometry of Li05. 



4.2.1 V & R-band photometry 

Comparing the pre-explosion photometry of the Li05 candidate progenitor with that of the coinci- 
dent source at late-time, we find that measurements in the V & R bands match within the errors. 

13 http://www.ast.cam.ac.uk/ wfcsur/technical/photom/colours/ 
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Figure 6. Archival pre-explosion and late-time post-explosion observations of the field of SN 2004et. (a&b): pre-explosion CFHT+CFH12K 
B&V observations, as studied previously by Li05; (c): pre-explosion CFHT+CFH12K R image of higher image quality and depth than previously 
studied; (d): pre-explosion INT+WFC sloan-i' image (e-h): WHT+AUX late-time BVRI observations taken ~ 3 years post-explosion when the 
SN has faded. The Li05 progenitor candidate is marked by the cross hairs. Note that this object is still visible post-explosion in BVR with similar 
magnitudes (see Table [HJ, proving that it is not solely the SN progenitor at these wavelengths. However, in the /-band the object is significantly 
brighter pre-than post-explosion. This excess is due to the unresolved progenitor star. All images are aligned and oriented such that North is up and 
East is to the left. 
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Table 7. Calibrated photometry of secondary standards in the field of SN 2004et. Photometry taken from lMisra et alj i2007l) and Li05. ID numbers 
correspond to those on Figure[8] Photometric errors are in parentheses. 



ID 


B 


V 


R 


/ 


1 


15.80(0.008) 


14.73(0.010) 


14.14(0.005) 


13.47(0.023) 


2 


15.51(0.008) 


14.15(0.009) 


13.41(0.004) 


12.63(0.023) 


3 


19.40(0.027) 


18.48(0.015) 


17.96(0.013) 


17.31(0.053) 


4 


19.44(0.040) 


18.29(0.030) 


17.60(0.020) 




5 


18.31(0.017) 


17.10(0.012) 


16.43(0.018) 


15.69(0.026) 


6 


15.56(0.008) 


14.28(0.009) 


13.63(0.004) 


12.81(0.023) 


7 


18.30(0.018) 


17.46(0.013) 


17.04(0.005) 


16.39(0.030) 


X 


18.43(0.018) 


17.66(0.013) 


17.28(0.043) 


16.67(0.040) 


9 


18.33(0.015) 


17.39(0.011) 


16.93(0.007) 


16.25(0.028) 



•* . 

— .* 



Figure 7. Pseudo-colour image created from WHT BVI observations of the site of SN 2004et, taken ~3 yrs post-explosion. The cross hairs mark 
the position of the object originally identified by Li05 as the candidate SN progenitor. This yellow object is still visible 3 years after the explosion. 
Oriented such that North is up and East is to the left. 

There are three possible solutions to these observations: 1) that the V & R flux of the pre- and 
post-explosion objects come from the same source or blended sources, with negligible contribu- 
tions from the progenitor and SN respectively; 2) the pre-SN V & R flux is entirely from the 
progenitor star while the late-time flux is from the SN; 3) the object consists of multiple blended 
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Figur e 8. Pre-explosion CFHT+CFH12K B image of the site of SN 2004et. Field standard stars with calibrated magnitudes taken from lMisra et alj 
12007 ) and Li05 are labelled (see Table[7}. North is up and East is to the left. 
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sources with the progenitor and SN contributing some significant fraction of the pre- and post- 
explosion V & R flux respectively. 

The second scenario seems the least likely as it requires the SN to be significantly brighter 



than expectec 



Maguire et al. 



(pred icted V ~26.1, and R ^26.2 by extrapolating the late-time light curves of 
201ol) while at the same time coincidentally matching the V&R magnitudes of its 



progenitor, and is ultimately ruled out by the HST and Gemini AO which resolves the object 
at the SN site into 3 distinct objects ( §14.31) . The third scenario still requires a certain amount 
of fine-tuning in order to match the pre- and post-SN photometry. However, the degree of fine- 
tuning can be reduced if we arbitrarily assume ever increasing flux contributions from the other 
constituent sources, so that it becomes possible to "wash-out" differences between the progenitor 
and SN magnitudes. The first scenario is the simplest, and assumes that the progenitor and the SN 
contribute negligible V&i?-band flux to the pre- and p ost-explosion image s respectively. However, 



Kotak et al. 



(|2009|) . and discussed in more 



late-time optical spectra of the SN site presented by 
detail in S14.3.11 clearly show broad emission lines that can only be attributed to the SN. Hence we 
know that SN 2004et does contribute at least some of the V&R flux in the late-time WHT images. 
These observations therefore point to the third scenario. 

4.2.2 B-band photometry 

In the 5-band observations it appears that the Li05 candidate may be brighter at late-time; B pre = 
24.02+0.21, B post = 23.62+0.06. If the difference is real we would assume that the SN contributes 
at least the excess post-explosion flux (B = 24.9+q?) with the remainder coming from neighbour- 
ing unresolved sources, or at most the entire post-explosion flux. We predicted the B magnitude 
of the SN at the time of the WH T observations ( ~ 3 y rs post-explosion) to be «24.7 from ex- 



trapolation of the SN light curve in [Maguire et al.l (|2010i) . This extrapolated magnitude is in good 



agreement with what we measure for the post-explosion excess, suggesting that 1) the SN is still 
visible, at least in B, and 2) the majority of the late-time flux comes from other unresolved sources 
that would also have been visible pre-explosion. These observations are in good agreement with 
what we infer from the V&J?-band photometry above. 

4.2.3 I -band photometry 

Photometry of the I-band observations shows that the pre-explosion source coincident with the 
Li05 candidate is significantly brighter than its late-time counterpart; I pre = 21.27 + 0.08, I post = 
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Table 8. Ground-based pre-explosion and late-time photometry of the source close to the position of SN 2004et 



Filter Pre-SN CFHT Pre-SN CFHT & INT Post-SN WHT 
(Li et al. 2005) (this work) 



B 24.30(0.48) 24.02(0.21) 23.62(0.06) 

V 22.93(0.31) 23.04(0.27) 23.06(0.04) 

R 22.50(0.18) 22.24(0.06) 22.27(0.05) 

I - 21.27(0.08) 21.99(0.06) 



21.99 + 0.06. Comparing this object to the surrounding stars, it is easy to see from Figure[6]that it 
is visibly brighter in the INT image than in the WHT frame. Such an excess in the pre-explosion 
image can logically be interpreted as a detection of the SN progenitor, which must have been 
intrinsically very red or heavily extinguished (or both) to have avoided detection in the BV&R- 
bands. An /-band magnitude for the unresolved progenitor can be trivially calculated from the 
difference between the pre- and post-explosion photometry; I prog = 22.06 + 0.20. 
We predicted the magnitude of the SN at the epoch of the WHT observations to be / «24.6, 



by extrapolating the /-band light curve of 



Maguire et al. 



(|2010r) . If correct, the predicted photom- 



etry would suggest that the SN contributes ~9 per cent of the WHT /-band flux, which should 
be discounted when calculating the progenitor magnitude. This would result in a slightly brighter 
progenitor magnitude of I prog = 21.96 + 0.20. However, in extrapolating the BVR&I SN pho- 
tometry we have assumed that there was no flattening or re-bright ening of these ligh t curves due 



Kotak et al. 



(|2009h found that, 



to light echoes or ejecta-CSM (circumstellar material) interaction, 
starting October 2007, SN 2004et did in fact re-brighten significantly in mid-IR Spitzer data. How- 
ever, previous mid-IR observations from August 2007 (aro und the time of our WHT observations) 
had still shown a declining light curve (|Kotak et al.ll2009|) . so our assumption above may not be 
unreasonable. We will discuss the IR-echo in m ore detail in §14.3.11 Later in this paper ( §14.4.21) 



Kotak et al 



( 2009 ) to constrain the true contri- 



we attempt to use the late-time optical spectra of 
bution of the SN to the late-time optical photometry, and subsequently to correct the progenitor 
magnitudes. 



4.3 High-resolution data: alignment and photometry 

Immediately apparent from the HST and Gemini AO images is that the previously unresolved 
source, coincident with the Li05 candidate progenitor, in fact consists of at least three distinct 
sources (see Figures l9T&fl0l). This provides conclusive evidence that the Li05 candidate identified 
in ground-based, seeing limited observations, was not solely the progenitor of SN 2004et. What 
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Figure 9. Each panel in this figure is centred on the SN position and oriented such that North is up and East is to the left, (a): Pre-CFHT if-band, 
(b&c): WFPC2 F606W and F814W and (d-f): NICMOS F110W, F160W and F205W late time observations taken almost 3 years post-explosion 
(2007 July 08). The single source close to the SN position in the ground-based image is resolved into at least 3 objects in the HST images (see 
Table [9] for photometry). One of these objects (indicated by the cross hairs) is exactly coincident with the SN position, but it is unclear if it is a 
detection of the SN alone. 



contribution the unresolved progenitor made to the pre-explosion images is the question we at- 
tempt to address in the remainder of §0] 

Photometry of the WFPC2 and NICMOS datasets was carried out using PSF fitting techniques 
within iraf daophot. Empirical PS F models were cr eated using several isolated stars in each image. 
The updated WFPC2 zeropoints of Dolphin! (|2000ar) were taken from Andrew Dolphin's websitd_, 
while NICMOS Vegamag zeropoints were determined using data from the STScI sita^. Flight 
system magnitudes measured for the three sources clustered close to the SN site are recorded in 
Tabled 

Alignment with the 2005 ACS/HRC F625W image of the SN confirmed that one of the three 
sources is exactly coincident with SN 2004et, having a measured separation of 3 + 8 mas. This 
object is labelled "Centre" in Table [9] and is indicated by the cross hairs in Figure |9j However, 
with the current dataset it is impossible to determine whether object Centre is solely a detection of 
the SN in all filters, or whether a second unresolved source (e.g. a close companion of the former 
progenitor star) is now also visible. 



14 http://purcell.as.arizona.edu/wfpc2_calib 

15 http://www.stsci.edu/hst/nicmos/performance/photometry 



The progenitors ofSNe 2006my, 2006ov & 2004et 29 



Table 9. WFPC2 and NICMOS late-time Vegamag photometry of resolved sources close to SN site. Each source is labelled with reference to its 
position relative to the SN position (see Figure |9j. Note the source labelled "Centre" is exactly coincident with the position of SN 2004et. 



WFPC2 NICMOS 
Object F606W F814W F110W F160W F205W 

Date 



East 

2007 Jul 08 23.88(0.04) 22.98(0.05) 

2008 Jan 19 23.90(0.03) 23.20(0.04) 
Centre 

2007 Jul 08 23.07(0.03) 22.80(0.04) 

2008 Jan 19 23.29(0.03) 23.06(0.04) 
South 

2007 Jul 08 >26.0 25.33(0.14) 

2008 Jan 19 >26.3 25.16(0.14) 



22.52(0.06) 21.67(0.05) 21.26(0.08) 

22.64(0.06) 21.95(0.06) >22.1 

22.17(0.06) 22.65(0.09) 21.70(0.15) 

22.42(0.06) 22.56(0.10) 21.22(0.09) 

22.86(0.07) 21.22(0.04) 20.78(0.05) 

22.89(0.07) 21.20(0.04) 20.75(0.05) 



What does appear clear is that object Centre is significantly brigh ter than one would ex pect the 



SN to be from simple extrapolation of the late-time light curves of LMaguire et al 



( 2010 ). At the 



time of the 2007 July 08 HST observations, we estimated the SN to have magnitudes of V «25.7 
and / «24.2. This compares with object Centre magnitudes of F606W = 23.07 + 0.03 (V-band = 
23.18) and FS14W = 22.8 + 0.03 (/-band = 2 2.78)B 



However, we know from 



Kotak et al. 



(|2009l) that the SN re-brightened significantly in the mid- 
IR beginning in October 2007. Although no mid-IR re-brightening was detected earlier than this 
date (including observations from August 2007) we cannot completely rule out the possibility of an 
optical light-echo, which could significantly increase the contribution of the SN to object Centre in 
the HST image from 2007 July 08. It is therefore also entirely possible that object Centre is solely 
a detection of SN 2004et. We con sider this case in §14.4. 11 while in 94.4.21 we use the late-time 



optical spectra of the SN site from 



Kotak et al 



(I2009() to constrain the minimum contribution of 



the SN to the late-time optical photometry. 



4.3. 1 NIR re-brightening ofSN 2004et and late-time optical spectra 

Although there is no direct evidence from the HST data for an optical light-echo associated with 
SN 2004et, we do observe a significant brightening of object Centre in the near-IR (NIR), both in 
the NICMOS F205W and Gemini Altair/NIRI i^T-band observations from 2008 Jan 19 and 2008 
Jul 10 respectively (see Figure ITOb. Object Centre, which is marginally detected in the F205W 
observation from 2007 July 08, increases in brightness by roughly 40 per cent by 2008 Jan 19. 
Photometry of the Gemini adaptive-optics observations was complicated by the large variation in 

16 HST photometry transformed to standard Johnson-Cousins V & I magnitudes using the colour transformations from 
http://purcell.as.arizona.edu/wfpc2_calib 
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Figure 10. Comparison of (a): NICMOS F205W image of the field of SN 2004et from 2007 July 08, (b): NICMOS F205W image from 2008 Jan 
08, and (c): Gemini Altair/NIRI K image from 2008 July 10. The images are centred on the SN position, which is also marked by cross hairs. The 
HST and Gemini images from 2008 clearly shows that SN 2004et has re-brightening significantly in F205W / K. 



the PSF across the FOV (due to anisoplanatic effects of the AO correction) and by having only one 
usable standard star (2MASS 20352132+6007102). Nevertheless, we estimated an approximate 
magnitude for object Centre in the 2008 July 10 Gemini image of K « 20.0. 

We attribute the increase in object Centre's F205W / .fiT-band flux to a re-brightening of SN 
2004et in the NIR, but stress once again that we cannot be certain that object Centre consists 
solely of the SN. Th e NIR re-brighteni ng of SN 2004et was contemporaneous with the mid-IR 
increase observed by lKotak et al.1 (|2009|) . 

Interestingly, we do not observe a brightening of object Centre in any of the other optical or 
NIR observations shortward of ~ 2yum. However, we note that the rate of decline in the optical 



bands between July 2007 and January 2008 is much slower than that measured by iMaguire et al 



(120101) for the earlier part of the radioactive tail; ~ 0.15 mag/100 days for the late-time HST pho- 
tometry compared to ~ 1.1 mag/100 days from Maguire et al. Furthermore, the late-time F160W 
photometry show no discernible decline whatsoever. All this could be interpreted as a flattening of 
the SN 2004et light curves caused by the same mechanism(s) that produced the F205W / K and 
mid-IR re-brightening. Another possible explanation is that a second unresolved source at the po- 
sition of object Centre (such as a close companion star to the former SN progenitor) is beginning 
to dominate over the fading SN. A combination of the two effects is also possible. Future HST 
observations of the SN site will be required to determine which is actually the case. 

Of immediate consequence for our study of the progenitor star in the optical, pre-explosion 
i mages are the late-ti me, optical spectroscopic observations of the site of SN 2004et presented 



in 



Kotak et al. 



200% . and reproduced here (Figure fTT)) with the permission of the author. The 
three spectra, taken with the Keck Telescope, span the period between ~2 to 3 yrs post-explosion, 
and clearly show broad emission lines that can only be attributed to the SN. The width of these 
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lines indicate expansion velocities of ~ 8500 km s~ l , and Kotak et al. interpret this as being due to 
the impact of the SN ejecta on the progenitor circumstellar medium. Such interaction would have 
resulted in a slower rate of decline (flattening) of the optical light curves. This lends support to the 
first scenario in the previous paragraph, but lacking further observations we still cannot rule out 
the possibility of a second unresolved source at the SN position. 

The Keck spectra are certainly contaminated by flux from other resolved sources close to the 
SN site. The spectra were reportedly taken under relatively poor conditions using a slit width of 
around 1.5" (R. Kotak, private communication). A glance at Figure [9] of this paper shows that, 
depending on the slit orientation, several objects could have contributed significant flux, not least 
the object immediately to the east of the SN position. It is therefore impossible to determine the 
total flux contribution made by SN 2004et to the late-time optical spectra / photometry. 



4.4 Progenitor luminosity and mass estimates 



In this section we attempt to derive luminosity and mass estimates for the progenitor of SN 2004et. 
In 34.2.3l we detected an unresolved progenitor star as an /-band flux excess in pre-explosion INT 
observations, whe n these were comp ared with WHT images taken ~3 yrs post-explosion. The 



optical spectra of 



Kotak et al. 



(|2009h (Fig. QJ]) show that the SN was still visible at this late- 
epoch, and hence we must correct for its contribution to the late-time optical photometry before 
calculating the magnitude of the unresolved progenitor star. 

As discussed above, determination of the total flux contribution of the SN to the late-time pho- 
tometry is impossible. The Keck spectra are contaminated by nearby, resolved sources, while ob- 
ject Centre in the HST observations may consist of more than just SN 2004et. We can, however, 
estimate the maximum and minimum contributions of the SN from the HST photometry and the 
Keck spectra. These two extremes are considered separately in the following subsections, where 
we also determine corresponding values for the progenitor photometry, luminosity and mass. 

We also note that there is some evidence from the photometry in Table [9] to suggest that object 
East may be variable. This adds an unknown systematic uncertainty to our calculations in the 
following subsections. Further HST observations of the SN site will be required to determine if 
object East is truly variable. 

To calculate the progenitor luminosity, and ultimately estimate its mass from stellar models, we 
require the distance to the host g alaxy, the extinc t ion alo ng the line of sight and the (approximate) 



metallicity of the progenitor star. 



Botticella et al. 



(|2009l) present a review of the distance estimates 
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to NGC 6946 in the literature, and employ the unweighted mean as their preferred value; /u = 
28.78+0.08 (d = 5.7+0.2 Mpc). We adopt the same value in this paper. 



The extinction of E(B - V) = 0.4 1 + 0.07 was estimated by 



Zwitter et al 



(I2004|) from the equiva- 



le nt width of the Na i D line s in a high resolution spectrum of SN 2004et, and using the calibration 



of 



Munari & Zwitten (| 19971) . The reddening uncertainty was chosen by Li05 to bracket the lower 



limit, which assumes no host-galaxy extinction. As pointed out by 



of extinction towards SN 2004et is estimated to be E(B-V) = 0.34 (ISchlegel et al. 



i05, the Galactic co mponent 



1998). 



We also attempted to constrai n the metallicity of the progenitor of SN 2004et from the available 



literature. 



Pilyugin et al. 



(|2004f ) studied 9 H n regions in NGC 6946, but none of these are close 
to the location of SN. We therefore used their calculated abundance gradient and the de-projected 
galactocentric radius of SN 2004et (Rg/Ris = 0.92) to estimate the metallicity at its position. This 
yielded an oxygen a bundance of 12+log(0/H) = 8.3 dex, which is consistent with that of the LMC 



( Hunter et al. 



2007). 



4.4. 1 Object Centre is SN 2004et 

Firstly, we consider the maximum possible flux contribution of the SN to the late-time photometry. 
This is the case where SN 2004et is assumed to contribute all the flux of object Centre in the late- 
time HST images (Figure[9]& Table©. In this scenario, the pre-explosion source (Li05 candidate) 
is logically assumed to consist of the other HST resolved objects (i.e. objects East and South) 
plus the SN progenitor. By subtracting the flux of objects East and South from the ground-based, 
pre-explosion photometry, we can calculate photometry for the progenitor star. To this end, the 
July 2007 F606W and FS14W photometry of object East was transformed to standard Johnson- 
Cousins V & I magnitudes using the colour transformations from Andrew Dolphin's website 11 . 
With no detection of object South in F606W, we were unable to transform its F&14W photometry 
to the standard system. However, we note that, in the case of the F814W-to-7 transformation, 
the correction is generally very small. Johnson-Cousins photometry for object East was therefore 
calculated as V = 24.25+0.03 and / = 22.94+0.04, while for object South we set / = FSUW = 
25.33+0.14. 

These magnitudes and their respective errors were converted to flux units, and subtracted from 
the V & /-band photometry of the pre-explosion source (Table [8]). This yielded progenitor star 
magnitudes of, V prog = 23.47+0.40 and I prog = 21.57+0.12. The apparent colour of the progeni- 
tor was therefore (V-I) = 1.90 + 0.41. Applying our best estimate of extinction (see 34.41) . we 
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calculated the progenitor's intrinsic colour; (V-I) a = 1.24 + 0.48. This colour range corresponds 
to supergiant stars of type F8 to K3, which are signifi cantly hotter than the M-t ype red super- 
giant one would expect for a type II-P SN progenitor. Drilling & Landoltl (|2000h show that the 
values of (V-I) + BC for such stars are quite similar, and we calculated an average of (V-I) 
+ BC = 0.75+0.15. Applying this to the /-band detection, along with values for the extinction 
and distance modulus detailed above ( 34.41) . we calculated the progenitor's bolometric magnitude, 
M bo i = -7.06 + 0.23 and its luminosity, log L/L Q = 4.72+0.09. By comparing this luminosity to the 
STARS stellar models of LMC metallicity in Figure [3] we estimated the progenitor initial mass to 
be 10+f Mo. 



4.4.2 Estimating the minimum contribution ofSN 2004et from late-time optical spectra 

In this section we have attempted to estimate the minimum flux contribution of the SN to the 
late-time photometry. In order to constrain the SN flux, we firstly utilised the sts das synphot 



Kotak et al. 



packag e to measure the SN emission line fluxes in the late-time optical spectra from 
(|2009|) . By interpolating the measurements from the 2007 April 13 and 2007 November 12 spectra 
we derived emission line magnitudes that were contemporaneous with our 2007 August 12 WHT 
photometry, assuming an exponential decay in flux. The reader should note that these ground-based 
spectroscopic observations were made under natural seeing conditions, and hence the continuum 
flux contains contributions from at least three different sources, that is those resolved in the HST 
and Gemini AO observations. 
Mag nitudes were first of all measured from the flux-calibrated spectra as presented in lKotak et al. 



(120091) (Figure dB, but the interpolated results, V spec = 22.22 and R spec = 21.59, were found to be 
significantly brighter than our WHT photometry (see Table [8]). V spec and R spec were brighter by 
factors of x2.2 and xl.9 respectively. Since we do not know the exact conditions under which the 
spectra were taken, such as seeing and slit orientation, we cannot easily quantify the contamina- 
tion due to other nearby resolved sources (see Figure [6]). In the R -band, at least, the blended source 
closest to the SN site is much brighter than any of the nearby objects. We suggest that the neigh- 
bouring resolved sources could not have resulted in the o bserved i?-band ex cess in the spectra, 
and subsequently conclude that the flux values reported by Kotak et al.1 (|2009|) are approximately 
a factor of 2 higher than their true values. In the analysis that follows we have re-scaled the spectra 
by a factor of x0.5, so as to be consistent with our WHT photometry. 



In order to estimate the the emission line magnitudes, the continuum of each spectrum was fitted 
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7000 
Wavelength (A) 



Figure 11. Late-time optical spectra of SN 2004et from lKotak et alj j2009h . taken at the Keck telescope. Fluxes are as presented by Kotak et al. 
The spectra have been scaled by the amounts shown in brackets. The upper two spectra have been shifted vertically for clarity. The fluxes are 
inconsistent with our WHT and HST late-time photometry and we estimate that the fluxes shown here are too larg e by approximately a factor of 2. 
at least for the 933d and 1 146d spectra. The broad emission lines, particularly H ff , are clearly visible. Courtesy of lKotak et 



oximately ; 
aI U2009i) 



and removed before the synphot task calcphot was run using the Landolt V & R passbands. Since 
calcphot resets all negative fluxes to zero, positive noise spikes lead to artificial brightening of 
the measured magnitudes. We therefore added a constant flux offset to the continuum subtracted 
spectra before running calcphot, which ensured that all fluxes were greater than zero. The offset 
magnitudes were measured separately, and then subtracted from the combined results to find the 
true emission line magnitudes. Interpolation of these results yielded estimates of the emission line 
magnitudes on the date of the WHT observation; 2007 August 12. Following this procedure we 
found that the SN emission lines could account for just 5 per cent of the total flux in V, but around 
25 per cent of that in the R -band. This is not surprising, since the R filter covers the wavelength 
region that includes the prominent H„ and [Ca II] emission features, while the V-band only grazes 
the Hq, profile. 



Kotak et al. 



(120091) were noisy and included 



In the /-band wavelength range the spectra from 
prominent sky lines. However, the [Ca II] emission feature lies within the bounds of this passband. 
We estimated its contribution to the WHT /-band observations in the same manner as described 
for the other filters above, except that in this case all flux values at wavelengths longer than 7600A 
were set to zero in order to suppress the sky lines. The [Ca II] emission was thereby estimated to 
contribute around 9 per cent of the total flux in /; a lower limit to the SN emission line flux in this 
filter. 

Of course, we expect the SN to have contributed flux to the continuum also, although the abso- 
lute level of this contribution is impossible to measure from the current dataset. By comparing the 
V & R magnitudes of the original 2007 April 13 and 2007 November 12 spectra, we found that 
around 14-18 per cent of the April flux is lost from the continuum by the time of the November 
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observation. These continuum losses equate to roughly 9x 10 4 mag d 1 . Between our WHT obser- 
vations on 2007 August 12, and the final spectrum on 2007 November 12, we calculate that the SN 
continuum would have faded by ~0.08 mag in V & R, again assuming an exponential decay in flux. 
This constitutes a lower limit for the continuum contribution of the SN in the WHT observations. 
Later spectra may well have revealed a further decline in the continuum due to the fading SN. 
However, lacking such extra information, we only consider flux that can be logically attributed to 
the SN using the current data, and reiterate that it is a minimum value. (The reader should note that 
the combined emission-line plus continuum contributions derived here are in excess of the decline 
observed between the 2007 July and 2008 January HST optical observations. Hence, no further 
constraint can be gleaned from the HST data). 

Our estimates of the minimum SN flux in V, R&I were subtracted from the WHT post-explosion 
photometry of the Li05 candidate in order to calculate the contribution from the other unresolved 
stars. Subtracting the result from the pre-explosion photometry of Li05 candidate left us with 
the flux of the progenitor star itself. In this way we estimated the magnitudes of the progenitor 
in the pre-explosion images to be V prog =25.\?> + ™ 2 , i? prog =23.39+0.26 and I prog =2l. 88+0.17. (We 
assume a similar continuum reduction in I as we have observed in V & R). Clearly, the calculated 
flux of V pr og is insignificant when compared to its error. Therefore, we use only R pr0 g and I prog to 
characterise the progenitor in this case. Using the extinction quoted in §14.4[ the intrinsic colour of 
the progenitor is calculated as (R-I) = 1.16+0.36. These colours correspond to supergiant stars of 
types K4 to M4. Again, we noted that the (V-I ) + BC values for thes e stars are quite similar, an d 



calculated an average value of 0.89+0.1 1 from LElias et al.1 (|1985l) and Drilling & Landoltl (|2000f) . 



The bolometric magnitude and luminosity were calculated from I prog to be M bol = -6.62 + 0.24, 
and logL/L Q = 4.54+0.10. This corresponds to an initial mass for the progenitor of 8+j M Q . 

Having corrected for the emission line fluxes and minimum continuum contributions of the SN, 
we conclude that these results are our best estimates for the spectral type, luminosity and initial 
mass of the progenitor of SN 2004et. However, we stress that future high-resolution observations 
of the SN site using HST will be required to determine if object Centre has faded further in order 
to confirm or adjust these results. 



5 CONCLUSIONS 

We have re-analysed the pre-explosion observations of the type II-P SNe 2006my, 2006ov and 
2004et. In the cases of SNe 2006my and 2006ov we have argued that the published progenitor can- 
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Table 10. Summary of photometry, luminosity and mass estimates for the progenitors of SNe 2006my, 2006ov and 2004et 



Supernova 


Photometry 


Distance Modulus 


Extinction 
E(B - V) 




Colour 


Spectral Type 


Luminosity 
logL/L 


Mass 
M 


2006my 


F814W>24.8 


31.74 ±0.25 


0.027 






MO I (assumed) 


<4.51 


< 13 


2006ov 


F814W > 24.2 


30.50 ± 0.40 


0.022 






M0 1 (assumed) 


<4.29 


< 10 


2004et (a) 
(b) 


7 = 21.57 ±0.12 
I = 21.88 ±0.17 


28.78 ± 0.08 
28.78 ± 0.08 


0.41 ± 0.07 
0.41 ± 0.07 


(V- 
(R- 


T)„ = 1.24 ±0.48 
/)„ = 1.16 ±0.36 


F8 1 to K3 1 
K4ItoM4I 


4.72 ± 0.09 
4.54 ±0.10 


10 + 5 
8* 



(a) assuming SN 2004et contributes all the flux of object Centre in late time HST images ( 34.4. U 

(b) SN 2004et contribution to object Centre estimated from late-time optical spectra ( jj4.4.2t 



didates (Li07) are most likely not coincident with their respective supernova sites in pre-explosion 
HST+WFPC2 observations. We concluded that the progenitors of SNe 2006my and 2006ov were 
probably not detected, and hence derived detection limits for each (see Table ITDT). Based on the 
assumption that the unseen proge nitors were red supergiants, a reasonable assumption given both 



SNe were classified as type II-P (|Chevalierlll976c lArnettill98Qh . we calculated upper luminosity 



and mass limits for the progenitors of SNe 2006my (log L/L Q < 4.51; m < 13M Q ) and 2006ov 
(log L/Lq < 4.29; m < 10M o ). Future high-resolution observations of each SN site will be required 
to confirm these results. 

Late-time WHT observations of the site of SN 2004et revealed that the yellow supergiant pro- 
genitor candidate, identified in pre-explosion CFHT BVR data by Li05, was still visible some 3 
years post-explosion. High-resolution, late-time Gemini AO and HST images resolve this object 
into at least 3 distinct sources, one of which (object Centre) is coincident with the SN. A hitherto 
unpublished, pre-explosion /'-band image from the INT was also analysed. Through comparison 
with the late-time WHT /-band data, we discovered a pre-explosion excess of flux coincident with 
the Li05 progenitor candidate. We attributed this to a detection of the unresolved progenitor star. 

Gemini Altair/NIRI and HST NICMOS observations revealed an apparent re-brightening of the 
SN in the NIR ~3.5 - 4 years post-explosion. This coincides with the re-brightening observed by 
Kotak et al. (2009) in the mid-IR. Late-time optical spectra presented by Kotak et al. confirmed 
that SN 2004et was still visible in the optical at the time of our WHT photometry. Correcting 
for the SN contribution we found an excess of pre-explosion flux in the /?-band, which we again 
attributed to the unresolved progenitor. By combining the R and /-band progenitor photometry 
we estimated the progenitor of SN 2004et to be a K4 to M4 supergiant, of logL/L Q = 4.54+0.10, 
and initial mass 8^ M (see Table ITOl). These constitute our best estimates for the spectral type, 
luminosity and mass of the progenitor of SN 2004et, but we stress that future observations may 
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show that the SN has faded still further, in which case our progenitor photometry would require 
revision. 
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APPENDIX A: A NOTE ON EXTINCTION 



In all cases in this paper, we have used our best estimate for the extinction towards each SN 
as a proxy for that of its progenitor star. We therefore assume that the extinction towards each 
SN is the same as that towards its progenitor before explosion. However, if the progenitor stars 
happ ened to be surrounded by a significant quantity of circumstellar dust, this might not be true 



(e.g., 



Botticella et al. 



explosion (|Dwek 



2009). Dust within a certain radius could be photo-evaporated by the SN 



1983|), effectively reducing the extinctio n measured towards a SN compared 



with that towards its progenitor star. 



Botticella et al. 



(|2009l) have estimated that it is possible for 



the UV-optical luminosity of a type II-P SN to destroy circumstellar dust equating to several tens 
of magnitudes of optical extinction. 



All this is of little concern if the progenitor star is free of s uch circumste 



Meikle et al. 



lar material, or if this 



20061) . If a progenitor 



material is far enough from the star to avoid evaporation (e.g., 
star was detected across a sufficient number and range of band-passes one could directly measure 
its extinction. Such direct measurement could reveal the presence of circumstellar dust that was 
subsequently destroyed by the SN. Unfortunately none of the three SNe studied in this paper fulfil 
these requirements. 

In §6.2 of lsmarttetalJ 



(|2009|) . we have discussed the progenitor extinction problem with re- 
lation to 20 type I I-P SNe, including the three featured in this paper. We refer the reader to 
Smartt et al. (2009J) for the detailed discussion, but highlight some of the main points here. The 



distribution of extinction values estimated for our progenitor sample w as more or less consi stent 



with that for optically selected red supergiants in the LMC and SMC (|Levesque et al 



2006J), al- 



beit with some evidence to suggest there were greater numbers of objects fro m our sample in th e 



Smartt et al 



2009) 



lowest extinction bin than would be typical for red supergiants (see Figure 5 
An arbitrary increase in A v of 0.3 mag for the 5 lowest reddening values in our sample brings 
the mean into line with that of the LMC/SMC red supergiants. The progenitors of SNe 2006my 
and 2006ov are two of the five objects ostensibly suffering the lowest extinctions. Applying the 
above arbitrary increase raises the upper mass limits for these two progenitors to: 2006my < 14M ; 
2006ov< 11M Q . 

Highly reddened (A v ~ 10) red supergiants do exist (e.g. Ivan Loon et alJl2005D a nd these, being 



optically obscured, would not be included in the sample of LLevesque et al.1 (|2006r) . Our compari- 
son sample is therefore biased towards lower extinctions, although to what extent is unclear since 
the relative number of heavily extinguished to low/moderately extinguished red supergiants is un- 
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Figure Al. The properties of the progenitor of SN 2004et, assuming arbitrary levels of extinction. Top panel: Range of extinction levels required 
to fit a star of certain spectral type to the progenitor photometry. Increased extinction can lead to higher luminosity/mass progenitors, but such 
stars are also of earlier/hotter spectral type. Bottom panel: The shaded area denotes the range of luminosities for a progenitor of a given spectral 
type/effective temperature. 

known. Without these relative numbers, we favour the progenitors to have come from a population 
of low/moderately extinguished red supergiants, although ultimately we cannot rule out that any 
of the three progenitors studied in this paper could have experienced much higher reddening than 
estimated. 

How would this essentially arbitrary source of extinction affect our results? Assuming arbitrary 
extinction, no constraints whatsoever can be placed on the mass limits or spectral types of the 
progenitors of SNe 2006my and 2006ov, since in both cases we do not detect a progenitor star. 
In the case of SN 2004et some constraints are possible. Assuming the progenitor photometry 
of 34.4.21 the colour of any fitted progenitor model must satisfy the observed (apparent) colour 
range of (R -/)= 1 .5 1 +0.31 . For a progenitor of given spectral type, a range of extinctions between 
minimum and maximum values will satisfy this condition. Hotter/bluer spectral types will require 
higher minimum extinctions, while such increases in extinction will rule out cooler spectral types. 
The result is that, as extinction is increased, we can fit progenitors of higher luminosity/mass, but 
these must also be of earlier spectral type. 

The results of a simple calculation using the photometry of the S N 2004et progenitor as der ived 
in 34.4.21 and the intrinsic (R—I) colours of supergiant stars from Drilling & Landoltl (|2000i) are 
plotted in Figure lATl The minimum and maximum values of extinction, A v , required to fit each 
spectral type are easily calculated from: 

A v = (A R -A 7 )/0.272 (Al) 
where, 

(A R -Aj) = (R-I)-(R-I) (A2) 
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The denominator of eqn . lAll comes from the ratio sA R /Ay =0.751 and Aj /A v =0 . 47 9 from 



Cardelli et al 



19890 . The minimum and maximum luminosities plotted in Figure ED were calculated from the 



I-band progenito r photometry (§14.4.21). the distance modulus ( §14.41) . the bolometric and colour 



corrections from 



Drilling & Landoltl (120001) . and the minimum and maximum extinction values 



derived for each spectral type, stars stellar models for stars of initial masses between 7-120 M 
are also plotted. These can be used as a rough guide to the initial stellar mass for a progenitor of 
given spectral type. Alternatively, the plotted luminosities can be compared directly with model 
endpoint luminosities in Figure [3] 



(Argo et al. 


2005; 


Misra et al. 


2007; 



Martf-Vidal et al 



2007b . while in §g3j] we show a NIR re- 



brightening some 4 years post-explosion and discuss the results of Kotak et al. (2009) who find 
evidence of ejecta-CSM interaction. These observations infer the presence of significant amounts 
of circumstellar material (CSM) around the progenitor star. One might therefore argue, with some 
justification, that the progenitor of SN 2004et was significantly extinguished by CSM, some of 
which may have been destroyed in the SN explosion. Our simple analysis suggests that if the 
extinction was greater than A v « 4.0, the progenitor star could not have been a red supergiant. On 
the other hand, if it was a red supergiant it could not have had an initial mass greater than ~16 M Q . 

It should be noted that to assume a progenitor of earlier spectral type automatically requires it 
to be more compact. For example a GO supergiant has a radius of R* « lOOR^. Such radii are 



not easy to reconcile with the fact that SN 2004et was a typical II-P (|Misra et al 



200 



2007 



Maguire et al 



Arnet 



3 l98cL whi 



Sahu et al. 



2010). SNe II-P require larger radii of the order 400- IOOORq (IChevalierlll97 



6 



which are typical for red supergiant stars. 



